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Relativistic Hartree-Fock-Bogoliubov (RHFB) theory with density-dependent meson-nucleon cou- 
plings is presented. The integro-differential RHFB equations are solved by expanding the different 
components of the quasi-particle spinors in the complete set of eigen-solutions of the Dirac equations 
with Woods-Saxon potentials. Using the finite-range Gogny force D1S as an effective interaction 
in the pairing channel, systematic RHFB calculations are performed for Sn isotopes and N = 82 
isotones. It is demonstrated that an appropriate description of both mean field and pairing effects 
can be obtained within RHFB theory with finite range Gogny pairing forces. Better systematics are 
also found in the regions from the stable to the neutron-rich side with the inclusion of Fock terms, 
especially in the presence of p-tensor couplings. 
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I. INTRODUCTION 

During the past decades, much success has been 
achieved in nuclear physics by relativistic density func- 
tional theories. One of the most outstanding candidates 
is the relativistic Hartrec approach with the no-sea ap- 
proximation, namely the relativistic mean field (RMF) 
theory 0, @, S, !, 11, S 0- Within the RMF frame- 
work, valuable information has been obtained for the 
structure of the nuclei in and far from the valleyof (3- 
stability, including both for ground states [1, 0,0] and 
excited states [6, 8]. At the same time, considerable ef- 
fort has been devoted to relativistic Hartree-Fock (RHF) 
theory @, EE EH El EE El EE El- However, because 
of its numerical complexity, for a long time, it failed in a 
quantitative description of nuclear systems. Only in the 
recent years, with the growth of computational facilities 
and the development of new methods, density dependent 
relativistic Hartree-Fock (DDRHF) theory has shown sig- 
nificant improvements in a quantitative description of nu- 
clear phenomena [H El, El, El El, E| with a similar 
accuracy as RMF. 

In DDRHF, the Lorentz covariant structure is kept in 
full rigor, which guarantees the self-consistent determina- 
tion of the spin-orbit coupling [ItJ and all well-conserved 
relativistic symmetries, e.g. , the pseudo-spin symmetry 
in the nuclear spectrum [18| . In addition, significant im- 
provements on the relativistic description of shell struc- 
tures have been gained with the newly introduced con- 
stituents by the Fock terms, i.e., the pseudo- vector pion 
and the p-tensor couplings. In Rcfs. 0, HE, the con- 
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sistency of the evolution of the shell structure has been 
considerably improved by the pion exchange potential, in 
fact, by its tensor part. With the inclusion of p-tensor 
couplings, the common disease of several artificial shell 
closures existing in the RMF calculations [24], has been 
cured in DDRHF theory and the pseudo-spin symmetry 
is also better preserved EE- Besides the Fock terms in 
the isovector channel, those derived from the isoscalar a 
and u) couplings are found to play a dominant role in re- 
producing the characteristic experimental Z-dcpcndcncc 
of the spin-orbit splitting around the sub-shell closure 
Z = 64 [IE- It has also been demonstrated that the 
isoscalar Fock terms are essential for self-consistent de- 
scription of the spin-isospin resonances within RPA [2lj | 
and prediction of neutron star properties [22l |. 

On the other hand, the development of the radioac- 
tive ion beam (RIB) facilities [26[ has opened a new 
frontier for nuclear physics, the field of exotic nuclei far 
from the valley of stability [H EL El El El, El, HE 
and the up grad es and constructions of the RIB facilities 
[34l , El, HE l37| in recent few years will provide us with 
new possibilities to study exotic modes in nuclear sys- 
tems. The current application of DDRHF is limited to 
nuclei in the /3-stability valley and pairing effects in open 
shell nuclei are treated only within the BCS approxima- 
tion E3, EE ■ I n weakly bound systems like exotic nuclei 
close to the drip lines, the Fermi surface of one type of 
nucleons is close to the particle continuum, and the single 
nucleon separation energies are comparable to the pair- 
ing gaps. This results in an enhancement of scattering of 
Cooper pairs into the continuum due to pairing correla- 
tions. Thus, it becomes necessary to include the contin- 
uous part of the single-particle spectrum to describe the 
unstable nuclei. 

It is now the general consensus that a unified and self- 
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consistent description of both mean field and pairing cor- 
relations can be obtained with the Bogoliubov transfor- 
mation and automatically the continuum effects are effi- 
ciently taken into account [38 , 39] . In this manuscript rel- 
ativistic Hartree-Fock-Bogoliubov (RHFB) theory with 
density dependent meson-nucleon couplings is presented 
as a natural extension of DDRHF. Its content is orga- 
nized as follows. In Section |TT] we introduce the gen- 
eral formalism of the RHFB theory with both zero-range 
(delta) and finite range (Gogny) pairing forces, where 
the integro-differential RHFB equations are solved by ex- 
panding the lower and upper components of the quasi- 
particle spinors on the complete set of solutions of the 
Dirac equation with a Woods-Saxon (DWS) type poten- 
tial. The comparison between different treatments of 
pairing correlations is discussed in Section IIIII and sys- 
tematic RHFB calculations are performed and discussed 
by taking Sn isotopes and N = 82 isotones as represen- 
tative cases. Finally, a brief summary is given in Section 

m 



II. GENERAL FORMALISM 

We briefly recall here the general features of DDRHF 
theory in order to make understandable its generalization 
to the RHFB case. More details can be found in Refs. [H], 
GjiGil whereas the effective interactions used in this work 
have been introduced in Refs. flTl [l9l. I20I . 



A. Energy Functional and Relativistic 
Hartree-Fock potentials 

As generally recognized, the nucleon-nucleon interac- 
tion is mediated by the exchange of mesons with isoscalar 
and isovector character. The understanding of nuclear 
structure at the microscopic level, therefore, has to be 
achieved in the same language. Consistent with this cri- 
terion, the model Lagrangian, i.e., the theoretical start- 
ing point, contains the degrees of freedom associated with 
the nucleon, the as u>-, p-, and 7r-meson fields, and the 
photon field (A) [l4j . Following the standard variational 
procedure of the Lagrangian [2, HH > one finds the equa- 
tions of motion for mesons, nucleons and photons, namely 
the Klein-Gordon, Dirac and Proca equations, and the 
continuity equation, i.e., the energy-momentum conser- 
vation relation, from which is derived the Hamiltonian 
of the system. In terms of the creation and annihilation 
operators (c^,c a ) defined by the stationary solutions of 
the Dirac equation, the Hamiltonian can be generally ex- 
pressed as, 



H 



a/3 aa'f3f3' 



c lc C 8 C 0' C a' VgBB'a'i C 1 ) 

4> 



where T a p represents the kinetic energy and the two- 
body terms V£*g, a , correspond to different types of me- 



son (or photon) nucleon couplings denoted by <j>, 

T al3 = [ dx$ a {x) (-17 . V + M) i>p (2) 



x D^(x,x')il)^{x')il) a -{x). (3) 

In the two-body interaction terms, the interaction ma- 
trices Y^{x,x') read as 



r a (x,x') = - g a (x)g a (x'), 
Tufax') = (g^) x idul^x' » 
T v p (x,x') = {g P i^f) x ■ (g P ^f) x , , 
<r i i\ _ 



(4a) 
(4b) 

(4c) 



r J(*>*') (f P ^ k rd k ) x ■ (f P v» l m) x , , (4d) 

T V^ X ') =^ (f P ° k ^d k ) x ■ (g P i,r) x , 

+ ^ I (9 P ^T) x -(f p a k »Td k ) x ,, (4e) 

T^x, x 1 ) = (Uf^d») x ■ (Ur l5lv d v ) x , , (4f) 

T A (x,x') = & - ( 7m (1 - t 3 )) x (y (1 - t 3 )) x , . (4g) 

In coordinate space, the propagators D < j,(x,x') for the 
meson fields have a Yukawa form 



D^x.x') 



1 e 



— mj, IX— X 



4tt \x — x' 



(5) 



For the photon field, the propagator D A (x,x') can be 
written as, 



D A (x,x') = 



1 1 

47T \X — X 1 



(6) 



In the above expressions (Eqs.(2-5)), M denotes the nu- 
cleon mass and m a (g a ), m w (g u ), m p {g p ,f p ), and 
(/tt) are the masses (coupling constants) corresponding 
to <7, cj, p, and tt mesons. In this paper, we use arrows 
to denote isospin vectors and bold types for vectors in 
coordinate space. 

In the Hamiltonian |T|), the indices a, (3, a', j3' run over 
all the single-particle states (ip a ) with positive energies 
(a = k) and negative energies (a = I). As it is commonly 
done in the mean field approach, the so-called no sea 
approximation is adopted and the contributions from the 
negative energy states are neglected. Then, the energy 
functional can be obtained from the following expectation 
value 



E = (<f> \H\<f> ), 



(7) 



where |$o) is the Hartree-Fock ground state in the no-sea 
approximation [l4| . In the energy functional ([7]) , the con- 
tributions of the two-body interactions consist of two 
parts, the direct (Hartree) and exchange (Fock) terms. 
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With only the direct contributions, Eq. (JT]) leads to the 
energy functional of the RMF theory. With both direct 
and exchange contributions we obtain the energy func- 
tional for the DDRHF theory. 

In spherically symmetric systems the Dirac spinor can 
be written as, 



1 ^G a (r)y l f ama {r) 

= - ; 

r \-Fa(r)y: rr Jr) 



(8) 



The radial wave functions G a {r) and F a (r) characterize 
the upper (large) and lower (small) components and y\ m 
are the spherical harmonic spinors. Here, the sub-index 
a = {a,m a } = {n a ,l a ,l' a , j a ,m a } contains the quantum 
numbers n a (number of nodes of the upper component 
Go), j ai m a (total angular momentum and its projection 
to the z-axis), and l a , l' a (orbital angular momenta with 
la + l' a — 2j'a). In the following, we will use the Latin 
indices for the sub-set {nll'j} and Greek indices for the 
full set {njll'm}. 

By taking the variation of the energy functional ([7]) 
with respect to the Dirac spinor ([8]) , we obtain the spher- 
ical Dirac Hartree-Fock equation as 



dr'h(r,r')ip(r') = eip(r) 



(9) 



where e is the single-particle energy (including the rest 
mass) and the single-particle Dirac Hamiltonian h(r,r') 
contains the kinetic energy ft, kln , the direct local potential 
hP and exchange non-local potential h E , 

/i kin (r, r') = [a-p + f3M] S(r - r'), (10a) 

h D (r, r') = [E T (r) 75 + E (r) + /3E s (r)] S(r - r'), 

(10b) 

' Y G (r,r>) Y F (r,r>) \ 
Jcfry) X F (r,r>) j ' 



/i E (r,r') 



(10c) 



In the above expression, the local self-energies Eg, 
Eo and E-r contain the contributions from the direct 
(Hartree) terms |, 1 i and the rearrangement 
terms [f|. The non-local self-energies Xg, X f , Yq and 
Yp come from the exchange (Fock) terms and they take 
the general form 

4tW) =E^S 2 (9^») r ^ b G (m^r,r') (g^G b ) r , , 

(11a) 

4tW) =E^t5'6 {9*Fb) r *Zr (m^ry) {g^F b ) r , , 

(lib) 

Ya?(rS) =E^t31 (g^M^im^ry) {g^G b ) r , , 

(11c) 

Y F t\r,r>) =E^at3'6 i3*G b ) r &%{mtfTy) ( 94> F b ) r , . 

(lid) 



In these expressions, g^ represents the coupling con- 
stants, % = y/2jl + 1, and ST^ h denotes the isospin fac- 
tors: 6 TaTb and 2 — S TaTb respectively for isoscalar and 
isovector channels. For example, one has M G = 3£ Xf = 
-M Yf = -M x ° = for the er-scalar coupling, and 

' 2 

^K^.r^^fe^) R LL (m a ;r,r'). (12) 



The prime on the sum in Eq. (|12|) indicates that L+l a +l b 
must be even, and Rl 1 l 2 stands for 



RL 1 L 2 (m l ;r,r') 



(z)K L2+h (z')6(z' - z) 



+ K LiH (z)I L2H (z')9(z - z')\ , (13) 

where z = m^r, and K L+ i are related to the spher- 

ical Bessel and Hankcl functions. The detailed expres- 
sions of all self-energies entering the HF potentials can 
be found in Ref. [Til ] excepted for the rearrangement po- 
tentials because the couplings there were assumed density 
independent. Here, the rearrangement potentials are of 
course included in the calculations. We observe that, in 
a non-relativistic reduction, the pion pseudo-vector cou- 
pling and the p-tensor coupling lead to central and tensor 
nucleon-nucleon interactions and therefore, they play a 
substantial role in determining the spin-orbit splittings 
and shell evolutions [H, [2(| ■ 

In realistic applications, one has to consider the nuclear 
medium effects. Within the RHF approach, some efforts 
have been devoted to considering the in-medium effects 
by introducing non-linear self couplings of the a and lo 



fields 15 1 or cubic and quadratic terms of the scalar field 



lq . Instead of the non-linear self couplings, here 
we assume a density dependence of the meson- nucleon 
couplings [10, HH E2] as we did before [l?], Gil, which 
looks more coincident with the model Lagrangian. 

As shown in Ref. [42j], the density dependence in 
mcson-nuclcon couplings leads to rearrangement terms 
E^, in the self-energy E^ in order to preserve the energy- 
momentum conservation, 



(14) 



For example, the rearrangement term due to the density 
dependence in cr-scalar coupling can be written as, 



4*) 



dga_ 
dp b 



g^r 2 



G b Y 



O) 



(<r) 



(15) 

where p s and p b are respectively the local scalar and bary- 
onic densities, and the Fock components X^ and Y b 
can be written as, 



Y, 



XI 



(<i>) 



= / dr' 



Y, 



G b 

(4>) y(« 



X G; X 



F b 



(r,r') 



(16) 
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B. Density-dependent Relativistic 
Hartree-Fock-Bogoliubov theory 

In open shell nuclei, the effects of pairing correlations, 
which lead to valence particles spreading over the or- 
bits around the Fermi level, have to be taken into ac- 
count, cither in the BCS approximation [43| or by the 
full Bogoliubov theory [44|. In terms of quasi-particles, 
the Bogoliubov theory unifies the treatment of ph- and 
pp-correlations in a self-consistent description of nuclear 
orbitals [45| . It is specially significant for the exploration 
in the regions far from the stability where the simple BCS 
method may break down. In the relativistic case 0, |47} 
earlier investigations within relativistic Hartree Bogoli- 
ubov (RHB) theory have shown that the scattering of the 
Cooper pairs into the continuum plays an important role 
for the formation of the neutron halos [H, [48| . Within 
the Bogoliubov scheme, the single-particle basis {c Q , cj,} 



and the quasi-particle basis {f3 a , } (a = 1, 
related by the following transformation 



= W 



pi 



ipv Vu 



Pa 

pi 



M) are 



(17) 



where ipu and tpy are the quasi-particle spinors, of the 
form of Eq. ([8]) in the spherical case. The transformation 
satisfies unitarity 



W f W = 1. 



(18) 



Following the standard procedure of the Bogoli- 
ubov transformation [44|, a relativistic Hartree-Fock- 
Bogoliubov equation can be derived as [46| . 



dr 



, /'/i(r,r / ) -A A(r,r') 
A(r,r') -h(r,r')+X 



Mr') 
ipy(r') 



E 



?Pu{r) 
ipv(r) 



(19) 



where the chemical potential A is introduced to pre- 
serve the particle number on the average. In the single- 
particle Hamiltonian h(r,r'), the retardation effects are 
neglected as is usually done in mean field calculations. 
The pairing potential can be written as, 



where the pairing tensor k is 

K a (r,r') = ip Va (r)*ip Ua {r'). 



(21) 



For the pairing interaction V pp in Eq. (|20[) . a phe- 
nomenological form is adopted as it has been done with 
great success in RHB theory 0, H3] and in conventional 
HFB theory 0, • The pairing force is either taken as 
a density-dependent two-body force in a zero range limit, 

V(r, r') - V S(r - r') \ (1 - a ■ a') (l - , (22) 

4 V Po J 



with an adjusted strength Vb, or as the pairing part of 
the Gogny force [5l[ , 

V(r,r') = V e(( r " r ')/^) 2 

i=i,2 (23) 

x (Wi + BiP° - H l P T - MiP a P T ) , 

with the parameters Hi, Wi, Bi, Hi and Mi (i = 1,2). 

In spherically symmetric systems the solution of the 
RHFB equations, i.e., the Dirac spinor ipu a and Vv Q can 
be written similarly to Eq. ((8|), 



, , if iG Ua {r)y l -: ma (r) 
r \-FvAr)y l 3 : ma {r) 
1 / iGv a (r)ytm a (r) 
r \-Fv a {r)y): m Sr) 



(24) 
(25) 



The RHFB equations (|19|) are then reduced to the system 
of coupled integro-differential equations, 



d 
dr 



r 



G Ua (r)-(E a + \--E_)F Ua (r) 



+ X Ua (r)+r J r'dr'A a {r,r')F Va {r') = 0, 

(26a) 

F Ua (r) + {E a + X-E+)G Ua (r) 



a n a 

2_, T 

dr r 



d 
dr 



r 



Y Ua {r)+r J r'dr'A a (r,r')G Va (r') = 0, 

(26b) 

G Va (r) + (E a -\ + Z-)F Va {r) 



+ X Va {r)+r J r'dr'A a {r,r')F Ua {r') = 0, 

(26c) 

F Va {r)-(E a -X + E + )G Va (r) 



a n a 

h T 

dr r 



- Y Va (r)+r J r'dr'A a (r,r')G Ua (r') = 0, 

(26d) 

where E a are the quasi-particle energies (without the rest 
mass), and the local self-energies £+ and £_ are 



T Zj S , 



£_ =£n - £<? - 2M. 



(27) 



In the radial RHFB equations (f2"o| , Xjj, Yu, Xy and Yy 
denote the contributions from the Fock terms, which are 
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of a general form similar to Eq. ([16 
dr 1 






(r,r') 



(r,r') 





(28a) 



(28b) 



For the nonlocal terms Xq. Xp, Yq and Yp above, one 
needs to replace the G and F components in Eqs. (TTTj) 
by the corresponding Gy and Fy in the general case, or 
by Gjj and F\j in the case of blocking. 

The pairing potentials A a (r,r') in Eqs. (f2l)|) can be 
expressed as 



A a (r,r') = -'£v: b p (r,r')K b (r,r'), 

b 

where the pairing tensor n(r, r') reads as 

H(r, r') J^i [G Ub {r)G Vb (/) + F Ub [r)F Vb {r')\ 



(29) 



1 



j 2 b [Gv b (r)G Ub (r') + F Vb (r)F Ub (r')} 



(30) 



Details of the pai ring interaction matrix element V£ b can 
be found in Ref. 13911. 



C. RHFB equations in Dirac Woods-Saxon basis 

In contrast to the RHB approach with <5-forces in the 
pairing channel where the radial equations (|26|) become 
differential equations, in RHFB theory the radial equa- 
tions are fully integro-differential. For zero-range (5-forces 
in the pairing channel the integral terms arise from the 
Fock terms, and for finite-range pairing forces they also 
come from the pairing channel. In coordinate space, it 
is difficult to solve such equations e.g., by a localization 
procedure similar to that adopted in Refs. 0, [TtJ- We 
therefore choose to solve them by an expansion of the 
Dirac-Bogoliubov spinors in an appropriate basis. 

In this work we solve the radial RHFB equations 
by using the Dirac Woods-Saxon (DWS) basis introduced 
by Zhou ct al. [52^ . This basis has been constructed 
for the investigation of weakly-bound nuclei. The set 
of DWS basis functions 



{[Eb,90(r,T)};e b ^ 0} , 



(31) 



are eigenfunctions (with eigenvalues £b) of a Dirac equa- 
tion with Woods-Saxon-likc potentials for Eo(r) ± T,s(r). 
They are determined by the shooting method in coordi- 
nate space within a spherical box of size i? max [EH . 



The U and V components of the Dirac Bogoliubov 
spinors (f24|) can be expanded as, 



N F N D 

i>u = ^2 U p9p + ^2 u d9d, 

p=l d=l 
N F No 

i'v = ^2 V p9p + X! Vd9d > 
p=l d=l 



(32a) 
(32b) 



where Np and Np> respectively correspond to the num- 
bers of positive (e p > 0) and negative (e^ < 0) energy 
states in the DWS basis. Obviously, because of spherical 
symmetry the quantum number k is preserved, i.e., the 
RHFB equations have to be solved for each value of k 
and the sums in the expansion (|32[) run only over states 
with the same n. For a fixed value of n we have the radial 
basis spinors 



9 P (r) = 



G P (r) 
F P (r) 



9d(r) 




(33) 



where the sub-indices p and d correspond to the number 
of nodes of the basis functions G p for positive energy and 
Fd for negative energy. 

In the DWS basis ((32]) the radial RHFB equations ([26]) 
are transformed to a matrix eigenvalue problem, 



(34) 



where H and A are (N F + N D ) x (N F + N D ) dimen- 
sional matrices, U and V are the column vectors with 
(Np + Np>) elements. From the expressions of the single- 
particle Hamiltonian h and pairing potential A given in 
the previous part we obtain the matrix elements of H 
and A as 




H 



Hi 



kin 



-7T + - I 

ar r 



drF n 



dr 



G n i , 



f dr[G n G n ,Z + + F n F n ,Z-] 

+ J dr(G n F n , + G n F n .) S T , 
dr I ' dr' (G n F„) 




(r,r') 




(35a) 



(35b) 



(35c) 



dr j dr'A K (r,r') 

x [G n {r)G n ,{r') + F n (r)F n ,(r')] , (35d) 



where n, n' run over the radial quantum numbers of the 
DWS basis states in Eq. (f3"3")) with both positive energies 
(n, n' = p) and negative energies (n, n' = d). 



G 



Before carrying out RHFB applications with the DWS 
basis, two constituents should be firstly decided, i.e., the 
size of the spherical box i? max and the number of states 
(Np and Nd) involved in the expansions ([52")) . In prac- 
tice, it is accurate enough to adopt the parameters of 
the DWS basis as R max = 20fm, N F = 28, N D = 12 
for the general applications whereas for weakly bound 
nuclear systems one needs to choose a larger spherical 
box radius (i? max == 24fm) and a larger number of states 
(N F =36). 



III. GENERAL APPLICATIONS OF THE RHFB 
THEORY 

We firstly examine the equivalence between different 
pairing mechanisms for stable nuclear systems. By using 
the parameter set PKA1 we perform the calcula- 
tions for the even-even Sn isotopes from 106 Sn to Sn 
by RHFB theory with Gogny and Delta pairing forces 
(referred to respectively by Gogny and Delta), an d by 
DDRHF with BCS pairing (denoted by BCS(5)) 
The comparisons are based on the fact that equivalent 
pairing gaps are obtained with different pairing treat- 
ments. For the DDRHF calculation with BCS pairing, it 
is performed completely in coordinate space [19( . 



TABLE I. Binding energy Eb/A and neutron radii r n for 
even-even Sn isotopes. The results are calculated by RHFB 
with Gogny and Delta pairing forces, and by DDRHF with 
BCS pairing [jg| |. in comparison with the experimental 
data [H3]. The used parameter set is PKA1 [l9( |. 







E/A (MeV) 






r„(fm) 




N 


Exp 


Gogny 


Delta 


BCS(<5) 


Gogny 


Delta BCS(<5) 


56 


-8.4327 


-8.4339 


-8.4423 


-8.4425 


4.456 


4.451 


4.470 


58 


-8.4688 


-8.4605 


-8.4687 


-8.4694 


4.508 


4.501 


4.523 


60 


-8.4961 


-8.4804 


-8.4877 


-8.4889 


4.558 


4.550 


4.573 


62 


-8.5137 


-8.4940 


-8.5000 


-8.5017 


4.606 


4.597 


4.620 


64 


-8.5226 


-8.5018 


-8.5063 


-8.5085 


4.651 


4.642 


4.665 


66 


-8.5231 


-8.5039 


-8.5071 


-8.5097 


4.695 


4.686 


4.708 


68 


-8.5166 


-8.5006 


-8.5029 


-8.5058 


4.735 


4.728 


4.748 


70 


-8.5045 


-8.4921 


-8.4937 


-8.4969 


4.772 


4.767 


4.785 


72 


-8.4879 


-8.4788 


-8.4799 


-8.4833 


4.805 


4.802 


4.818 


74 


-8.4674 


-8.4613 


-8.4616 


-8.4652 


4.835 


4.834 


4.847 


76 


-8.4436 


-8.4401 


-8.4396 


-8.4431 


4.863 


4.863 


4.874 


78 


-8.4168 


-8.4157 


-8.4145 


-8.4175 


4.889 


4.889 


4.897 


80 


-8.3869 


-8.3882 


-8.3871 


-8.3889 


4.913 


4.913 


4.917 


82 


-8.3549 


-8.3579 


-8.3579 


-8.3579 


4.935 


4.935 


4.935 


84 


-8.2779 


-8.2752 


-8.2744 


-8.2733 


4.993 


4.991 


5.001 


86 


-8.1990 


-8.1934 


-8.1916 


-8.1900 


5.050 


5.046 


5.062 



In Table [I] are shown the binding energy Eb/A and 
neutron radii r„, extracted from the calculations with 
Bogoliubov and BCS pairings. From TableUone can find 
good agreement on the binding energies since the studied 



nuclei are located in the stability valley. For the neutron 
radii, there exist some minor systematic deviations be- 
tween the results of Bogoliubov and BCS pairings. Ex- 
cept for the magic nuclei, the calculations with BCS pair- 
ing present slightly larger values (~ O.Olfm) than those 
given by Bogoliubov pairings. 
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FIG. 1. Neutron canonical single-particle energies for Sn, 
calculated by RHFB with Gogny and Delta pairing forces, 
and by DDRHF with BCS pairing [T^j. Horizontal error bars 
denote the occupation probabilities of the states and filled 
circles represent the Fermi levels. See the text for details. 

Taking 124 Sn as an example, in Fig. [T] arc shown the 
neutron canonical single-particle energies and the occu- 
pation probabilities (in horizontal error bars) extracted 
from RHFB calculations with Gogny and Delta pair- 
ing forces. For comparison are also shown the results 
from DDRHF calculations with BCS pairing. In the 
Bogoliubov scheme the canonical single-particle states, 
i.e., the eigenstates of the density matrix, can be ob- 
tained with the canonical transformation from the Bo- 
goliubov quasi-particle to the canonical basis (4f|. With 
the BCS approximation the density matrix and single- 
particle Hamiltonian do commute. The corresponding 
single-particle energies are therefore the canonical ones. 

As shown in Fig. [1] there is no distinct difference in 
the occupation probabilities (denoted by horizontal error 
bars) between different pairing treatments because of the 
existence of the shell gap 82. For the single-particle ener- 
gies, the calculations with Bogoliubov and BCS pairings 
provide identical values for the states below the Fermi 
level. For the states above, particulary the low-/ ones, 
remarkable deviations are found. As seen from the oc- 
cupation densities in Fig. [21 different pairing treatments 
lead to identical radial distributions for the deeply bound 
v2p states. For the v3p states, the occupation densities 
given by BCS calculations become rather diffuse at large 
distance although they are weakly bound. In contrast the 
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FIG. 2. The occupation densities of the v2p and v3p states, extracted from RHFB calculation with Gogny and Delta pairing 
forces, and from DDRHF (BCS(<5)) with BCS pairing [t<J. See the text for details. 



calculations with Bogoliubov pairings still present appro- 
priate asymptotic behavior at large distance even when 
the states lie beyond the particle continuum threshold. 
From Fig. [5] one may recognize that within the Bogoli- 
ubov scheme the occupation densities are properly local- 
ized inside the nucleus such that the continuum effects 
can be efficiently taken into account. For the stable nu- 
clei, this is less important, e.g., in 124 Sn the scattering of 
Cooper pairs into the continuum is blocked by the shell 
gap 82. In the weakly bound nuclei the valence orbits 
may gather around the particle continuum threshold and 
the continuum effects are then strongly enhanced. As 
shown in Fig. [2j such effects can be self-consistently and 
efficiently taken into account by the Bogoliubov transfor- 
mation jH, [3^] . 

In the applications of the RHFB theory with the zero- 
range pairing force, the cut-off on the quasi-particle en- 
ergy is an important ingredient as well as the pair- 
ing strength Vq. In the above calculations the pairing 
strength is set to Vq = 325MeV with the quasi-particle 
energy cut-off ~ lQOMeV. Compared to the zero-range 
pairing force, the finite range Gogny force is of less ar- 
bitrariness because of the finite range and natural cut- 
off. In addition, an appropriate description of the mean 
field can also be provided by the Gogny force in the non- 
relativistic calculations and therefore better systematics 
is expected with the Gogny- type pairing force. 

Now we aim for the systematical study of both pairing 
correlations and mean fields by considering Sn isotopes 



from 100 Sn to 137 Sn, and N = 82 isotones from 129 Ag to 
153 Lu as representatives. The calculations use the RHFB 
theory with the parameter sets PKA1 [l^| (with p-tensor 
couplings) and PKOl [TtJ (without p-tensor couplings), 
and they are compared to those obtained by RHB theory 
with the parameter set DD-ME2 (55|, one of the most 
successful candidates in the existing RMF effective inter- 
actions. In the following, the finite range Gogny force 
D1S [5l| is adopted in pairing channel. For the iso- 
topes (isotones) with an odd neutron (proton) number, 
the blocking effects have to be taken into account. In 
the corresponding calculations, we blocked different or- 
bits around the Fermi surface, which can be provided by 
the calculations of the neighboring even isotopes or iso- 
tones, and we chose the state with the largest binding 
energy \Eb\ as the ground state. 

In Table HT1 and Table HTT1 we show the binding energies 
per particle Eb /A for Sn isotopes and N = 82 isotones, 
respectively, as well as the blocked orbits (jb) for the 
odd-A isotopes. For the odd Sn isotopes we find, except 
for 123 Sii73, the same blocking configurations for the pa- 
rameter sets PKOl and DD-ME2, which provide similar 
neutron spectra, e.g., for 132 Sn (see Ref. [lj|). However, 
PKA1 shows very different blocking results for N < 65. 
This is mainly due to the fact that the pseudo-spin part- 
ners (^l<77/2, ^2g? 5 / 2 ) near the Fermi surface are somehow 
degenerate in the results of PKAl [l9[. In contrast the 
calculations with PKOl and DD-ME2 present remark- 
able gaps between these two states, i.e., the artificial 



8 



TABLE II. The binding energies per particle Eb /A (MeV) of Sn isotopes and the blocked neutron (u) orbits jb of the odd 
isotopes. The results are calculated by RHFB with PKA1 [!§] and PK01 [JJ], RHB with DD-ME2 [El, in comparison to the 
data The quantities A are the r.m.s. deviations. 





Exp 


PKA1 


PKOl 


DD-ME2 




Exp 


PKA1 


PKOl 


DD-ME2 


N 


E B 1 A 


E B 1 A 


E B /A 


E B /A 


N 


E B /A 


&Bj A 


jb 


&B/A 


jb 


fT 1 A 

EjBjA 


jb 


50 


— 8.2479 


—8.3097 


—8.2831 


—8 2635 


51 


—8.2740 


—8.3242 


u2d 5 / 2 


—8.3027 


flg 7 / 2 


—8.2804 


ulg 7/2 


52 


—8.3244 


—8.3587 


—8.3454 


—8.3198 


53 


—8.3420 


—8.3686 


u2d 5 / 2 


—8.3588 


vlg 7 /2 


—8.3327 


ulg 7 /2 


54 


—8 3836 

lj. 


—8.4001 


—8 3969 


—8 3688 


55 


—8 3965 


—8.4047 


v2d 5 / 2 


—8.4046 


vlg 7/2 


—8.3778 


vlg 7 / 2 


56 


—8.4327 


—8.4340 


—8.4390 


—8.4109 


57 


—8.4401 


—8.4327 


v2d 5/2 


—8.4413 


v2d 5/2 


—8.4158 


v2d 5/2 


58 


—8.4688 


—8.4606 


—8.4724 


—8.4463 


59 


—8.4706 


—8.4551 


vlg 7 /2 


—8.4715 


v2d 5 / 2 


O A A Ci l~7 

—8.4487 


v2d 5 / 2 


60 


—8.4961 


—8.4805 


—8.4977 


—8.4746 


61 


—8.4932 


O A TO O 

—8.4733 


v3Sx/2 


—8.4928 


v2d 5/2 


—8.4740 


v2d 5/2 


62 


—8.5137 


—8.4942 


—8.5149 


—8.4957 


63 


—8 5069 

L' ■ J V / V J t / 


—8.4854 


v3Sx/2 


—8.5049 


v2d s/2 


—8.4898 


v2d 5 / 2 


64 


—8 5226 


—8.5019 


—8.5243 


—8 5085 


65 


—8.5141 


—8.4912 


v3Sx/2 


—8.5117 


u3s 1/2 


O A C\C\T 

—8.4997 


u3s 1/2 


66 


—8.5231 


—8.5041 


—8 5260 


—8.5122 


67 


—8 5096 


—8.4907 


v3sx/2 


—8.5117 


u3a 1/2 


—8.5010 


u3s 1/2 


68 


—8.5166 


—8.5007 


—8.5213 


—8 5080 


69 


—8.4995 


O A O A A 

—8.4844 


v3Sx/2 


—8.5046 


u3s 1/2 


—8.4935 


u3s 1/2 


70 


—8.5045 


—8.4922 


—8.5110 


—8.4976 


71 


—8.4853 


—8.4725 


v3Sx/2 


—8.4920 


u3s 1/2 


—8.4793 


u3s 1/2 


72 


—8.4879 


—8.4790 


—8.4960 


—8.4820 


73 


—8.4673 


-8.4574 


vlhn/2 


-8.4744 


u3s 1 / 2 


-8.4610 


Vlhn/2 


74 


-8.4674 


-8.4615 


-8.4768 


-8.4624 


75 


-8.4456 


-8.4391 


Vlhxl/2 


-8.4529 


vlh ll/2 


-8.4406 


vXhxx/2 


76 


-8.4436 


-8.4403 


-8.4536 


-8.4395 


77 


-8.4208 


-8.4171 


vlh 11/2 


-8.4285 


vlhu/ 2 


-8.4169 


ulh 11/2 


78 


-8.4168 


-8.4158 


-8.4265 


-8.4139 


79 


-8.3928 


-8.3917 


vlh 11/2 


-8.4002 


Vlhn/2 


-8.3905 


ulhu/2 


80 


-8.3869 


-8.3883 


-8.3956 


-8.3858 


81 


-8.3629 


-8.3633 


vlh-n/2 


-8.3677 


Vlhn/2 


-8.3618 


vXhxx/2 


82 


-8.3549 


-8.3580 


-8.3605 


-8.3556 


83 


-8.3107 


-8.3103 


v2f 7/2 


-8.3093 


^2/7/2 


-8.3034 


v2ft/ 2 


84 


-8.2779 


-8.2754 


-8.2757 


-8.2644 


85 


-8.2320 


-8.2277 


v2f 7/2 


-8.2246 


^2/7/2 


-8.2123 


v2f 7/2 


86 


-8.1990 


-8.1936 


-8.1921 


-8.1744 


87 


-8.1530 


-8.1455 


u2f 7/2 


-8.1413 


v2f 7/2 


-8.1222 


v2f 7/2 


A 




0.0197 


0.0115 


0.0137 


A 




0.0177 


0.0095 


0.0146 



TABLE III. Same as Table HU but for N = 82 isotones 





Exp 


PKA1 


PKOl 


DD-ME2 




Exp 


PKA1 


PKOl 


DD-ME2 




Eb/A 


Eb/A 


Eb/A 


Eb/A 




Eb/A 


Eb/A 


jb 


Eb/A 


jb 


Eb/A 


jb 


130 Cd 


-8.2561 


-8.2563 


-8.2699 


-8.2491 


129 Ag 


-8.1930 


-8.1887 


71T39/2 


-8.2046 


7lTfl9/2 


-8.1772 


1T39/2 


132 Sn 


-8.3549 


-8.3580 


-8.3605 


-8.3556 


131 In 


-8.2988 


-8.2989 


7T139/2 


-8.3056 


7lTfl9/2 


-8.2941 


71T39/2 


i34 Te 


-8.3838 


-8.3818 


-8.3998 


-8.3888 


133 Sb 


-8.3649 


-8.3626 


nlg 7/2 


-8.3730 


TlflT/2 


-8.3659 


1T37/2 


136 Xe 


-8.3962 


-8.3912 


-8.4208 


-8.4063 


135j 


-8.3848 


-8.3788 


nlg 7/2 


-8.4031 


lT.97/2 


-8.3911 


71T37/2 


138 Ba 


-8.3934 


-8.3869 


-8.4253 


-8.4089 


137 Cs 


-8.3890 


-8.3807 


TTlgr/2 


-8.4157 


Ttlgr/2 


-8.4009 


71T37/2 


140 Ce 


-8.3764 


-8.3694 


-8.4123 


-8.3956 


139 La 


-8.3781 


-8.3685 


nlg 7/2 


-8.4107 


■Klg 7/2 


-8.3951 


7rlg 7/2 


142 Nd 


-8.3461 


-8.3395 


-8.3787 


-8.3618 


141p r 


-8.3540 


-8.3453 


n2d 5/2 


-8.3869 


n2d 5 / 2 


-8.3715 


Tv2d 5/2 


144 Sm 


-8.3037 


-8.2979 


-8.3312 


-8.3140 


143 Pm 


-8.3178 


-8.3097 


Tl2d 5 /2 


-8.3464 


n2d 5 / 2 


-8.3305 


n2d 5 /2 


146 Gd 


-8.2496 


-8.2449 


-8.2723 


-8.2548 


145 Eu 


-8.2693 


-8.2613 


7r2d 5 / 2 


-8.2922 


7t2c?5/ 2 


-8.2759 


n2d 5 /2 


148 Dy 


-8.1809 


-8.1810 


-8.2032 


-8.1853 


147 Tb 


-8.2067 


-8.2022 


7r2d 3/2 


-8.2268 


7r2d 5/2 


-8.2100 


Tv2d 5/2 


i5o Er 


-8.1022 


-8.1074 


-8.1250 


-8.1065 


149 Ho 


-8.1335 


-8.1346 


7t3si/ 2 


-8.1528 


■Klhxl/2 


-8.1356 


rrlftn/2 


152 Yb 


-8.0157 


-8.0252 


-8.0384 


-8.0196 


151 Tm 


-8.0501 


-8.0575 


7t3si/ 2 


-8.0710 


nlhxx/2 


-8.0533 


nlhxx/2 


154 Hf 


-7.9180 


-7.9354 


-7.9442 


-7.9250 


153 Lu 


-7.9593 


-7.9719 


7t3Si/2 


-7.9810 


■Klhxx/2 


-7.9629 


nlhxx/2 


A 




0.0071 


0.0247 


0.0099 


A 




0.0071 


0.0222 


0.0100 



shell closures N = 58 [HHH. In Table [U a long-range 
blocking is found in vs\i-z (more than 4 odd isotopes), 
which implies that the low-^ states are more favored by 
the blocking effects. For the odd N = 82 isotones we 
find in Table IIIII identical blocking on the neutron rich 
side (Z < 63: 145 Eu) for PKAl, PKOl and DD-ME2. 



When Z > 65 ( 147 Tb), PKAl gives a different blocking, 
e.g., the blocking favored state nlsi/2- In the last rows 
of Table [Til an d Table IIIII we show the the root mean 
square deviations A (averaged over the isotopes in the 
column) of the binding energy Eg /A from the experimen- 
tal values [iH for both even and odd nuclei. They indi- 
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52 56 60 64 68 72 76 80 84 88 48 52 56 60 64 68 72 

N Z 



FIG. 3. (color online) Single-neutron separation energies (MeV) along Sn isotopic (S n : left panels) and N = 82 isotonic (S p : 
right panels) chains. The results are calculated by RHFB with PKA1 [3, PKOl 0, and by RHB with DD-ME2 [H, in 
comparison to the experimental data [54| . 



N n 

52 56 60 64 68 72 76 80 84 88 52 56 60 64 68 72 76 80 84 88 




FIG. 4. (color online) Two-nucleon separation energies (MeV) (left panels) and the deviations (right panels) from the exper- 
imental data [Hj for Sn isotopes (Sa^: upper panel) and N — 82 isotones (S2 P '- lower panel). The results are calculated by 
RHFB with PKA1 [l§] and PKOl 0, and by RHB with DD-ME2 [H]. See the text for details. 



cate that the three models, RHF with p-tensor couplings 
(PKA1), RHF without p-tensor couplings (PKOl), and 
RMF (DD-ME2), present comparable quantitative accu- 
racies, and PKOl provides the best overall agreement for 
the Sn isotopes whereas PKA1 presents the best overall 
descriptions for N = 82 isotones. 

From the binding energies in Tablc|TT]and Table HTT1 we 
have extracted the singlc-nucleon and two-nucleon sep- 
aration energies to study the systematics of both mean 



fields and pairing correlations. Fig. [3]presents the single- 
neutron separation energies S n of Sn isotopes from 101 Sn 
to 138 Sn (left panels) and the single-proton separation en- 
ergies S p of N = 82 isotones from 130 Cd to 153 Lu (right 
panels), in comparison with the experimental data from 
Ref. [5J]. It is well known that the odd-even differences 
on the single-nucleon separation energies reflect the ef- 
fects of the pairing correlations. In Fig. [31 PKA1, PKOl 
and DD-ME2 present comparable and satisfactory quan- 
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titative agreements with the data for both isotopic and 
isotonic chains, which means that the appropriate de- 
scription of the pairing correlations can be provided by 
the RHFB theory with the finite-range Gogny pairing 
force. From Fig. [3l one can find some systematics in 
the results of these three models. On the neutron-rich 
side, i.e., after 132 Sn for Sn isotopes and before Gd for 
N — 82 isotones, PKA1 shows a better agreement than 
PKOl and DD-ME2. On the proton rich side, these three 
models show similar accuracy. 

In Fig. 0] are shown the two-nucleon separation en- 
ergies (plot a and b) and the deviations (plot c and d) 
from the experimental data for Sn isotopes (plot a and c) 
and N = 82 isotones (plot b and d). It can be seen that 
PKA1, PKOl and DD-ME2 reproduce well the data in a 
rather wide range, the deviations being within ±0.5 McV. 
As we know, the sudden change on the two-nucleon sep- 
aration energy in general reflects the existence of signifi- 
cant structure (e.g., at 132 Sn). From Fig. 0J; and Fig. 2H, 
one can see that PKA1 shows a different agreement from 
PKOl and DD-ME2. Along the Sn isotopic chain, PKAl 
presents good quantitative agreement from N = 61 to 87 
and large deviations are found on the proton rich side. In 
the results calculated by RHB with DD-ME2, large de- 
viations are seen on both neutron and proton rich sides 
as shown in Fig. Hfc- Among these three effective inter- 
actions, PKOl provides the best overall agreement with 
the data for Sn isotopes while for N = 82 isotones (Fig. 
@Ji) PKAl presents the best overall agreement. 

Concerning the separation energies, better systemat- 
ics are obtained from the stable region to the neutron 
rich side with the inclusion of Fock terms, especially with 
the presence of p-tensor couplings, e.g., around 132 Sn in 
Sn isotopic chain as well as the region around 140 Ce in 
N = 82 isotones (see right panels of Fig. g]). In fact, 
such improvements are consistent with the elimination of 
the artificial shell closures 58 and 92 0, [24| beyond the 
magic gaps 50 and 82, which may change the mean fields 
and pairing effects. These artificial shell closures appear 
in all RMF models, and in RHF they can be eliminated 
with the inclusion of the p-tensor couplings. In addition, 
the improved systematics from the stable region to neu- 
tron rich side are meaningful for the reliable exploration 
of the nuclear systems with extreme neutron-to-proton 
ratios. 

In principle, with the model Lagrangian based on 
meson-exchange nucleon-nucleon interactions one could 
have the same degrees of freedom in RMF as in RHF. 
However, the pion pseudo-vector and rho-tensor cou- 
plings cannot be efficiently taken into account by the 
RMF because of the lack of exchange terms. As pointed 
out in Refs. 0, [2(J, these two couplings bring indeed 



significant improvements on the description of the shell 
structure and its evolution while because of their nature, 
they do not bring much additional freedom to the de- 
scription of binding energies. This is the reason why three 
different models provide equivalent accuracy on the bind- 
ing energies of Sn isotopes and N = 82 isotones. Even 
though, distinct deviations still exist between RHF and 
RMF, or between RHFB and RHB, in the systematic 
behaviors of the binding energies. 



IV. SUMMARY 

In this paper, we have introduced the relativistic 
Hartree-Fock-Bogoliubov (RHFB) theory with density- 
dependent meson- nucleon couplings. The RHFB equa- 
tions are solved by an expansion of the Dirac-Bogoliubov 
spinors on a relativistic Dirac Woods-Saxon (DWS) basis. 
By taking the finite range Gogny force D1S as the pairing 
force, we have performed RHFB-DWS calculations for 
both stable and weakly bound nuclei. The parameters of 
the DWS basis are determined for the applications of the 
RHFB theory in exotic as well as stable nuclei. The quan- 
titative agreement between Bogoliubov and BCS pairings 
in describing the stable open shell nuclei was shown by 
taking the even Sn isotopes as the representatives. We 
have applied the RHFB theory with the Gogny pairing 
force to the study of Sn isotopes and N = 82 isotones, 
and demonstrated that the RHFB theory with the finite- 
range Gogny force in the pairing channel can provide an 
appropriate quantitative description of both mean field 
and pairing correlation effects. In addition, better sys- 
tematics from the stable region to the neutron-rich side 
are obtained with the inclusion of Fock terms, especially 
with the presence of p-tensor couplings which can elim- 
inate artificial shell closures at 58 and 92. In fact, such 
improvements on systematics are meaningful for reliable 
explorations of exotic regions. 
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